Optomechanical scheme for the detection of weak impulsive forces 

David Vitali^, Stefano Mancini^^^, and Paolo Tombesi^ 
^ INFM, Dipartimento di Matematica e Fisica, Universita di Camerino, 1-62032 Camerino, Italy 
^INFM, Dipartimento di Fisica, Universita di Milano, Via Celoria 16, 1-20133 Milano, Italy 

(Received: February 1, 2008) 

We show that a cooling scheme and an appropriate quantum nonstationary strategy can be used 
to improve the signal to noise ratio for the optomechanical detection of weak impulsive forces. 
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A mechanical oscillator coupled to an optical mode by 
the radiation pressure provides a sensitive device able 
to detect very weak forces. Relevant examples are in- 
terferometers for the detection of gravitational waves |^] 
and atomic force microscopes [Q. Up to now, the ma- 
jor limitation to the implementation of sensitive optical 
measurements is given by thermal noise |3| . It has been 
proposed in Ref. Q to reduce thermal noise by means of 
a feedback loop based on homodyning the light reflected 
by the oscillator, playing the role of a cavity mirror. The 
proposed scheme is a sort of continuous version of the 
stochastic cooling technique used in accelerators (|], be- 
cause the homodyne measurement provides a continuous 
monitoring of the oscillator's position, and the feedback 
continuously "kicks" the mirror in order to put it in its 
equilibrium position. This proposal has been experimen- 
tally realized in Ref. using the "cold damping" tech- 
nique raL which is physically analogous to that proposed 
in Ref. and which amounts to applying a viscous feed- 
back force to the oscillating mirror. 

Both the "stochastic cooling" scheme of Ref. Q and 
the cold damping scheme of Ref. ^ cool the mirror by 
overdamping it, thereby strongly decreasing its mechan- 
ical susceptibility at resonance. As a consequence, the 
oscillator does not resonantly respond to the thermal 
noise, yielding in this way an almost complete suppres- 
sion of the resonance peak in the noise power spectrum, 
which is equivalent to cooling. However, the two feed- 
back schemes cannot be directly applied to improve the 
detection of weak forces. In fact the strong reduction of 
the mechanical susceptibility at resonance means that the 
mirror does not respond not only to the noise but also to 
the signal, and we shall see that the signal to noise ratio 
(SNR) of the device in stationary conditions is actually 
never improved. Despite that, here we show how it is 
possible to design a nonstationary strategy able to sig- 
nificantly increase the SNR for the detection of impulsive 
classical forces acting on the oscillator. This may be of 
crucial importance in the field of metrology [p| , as well 
as for the detection of gravitational waves |^ We use 
a quantum treatment, allowing us to show why a clas- 
sical approach provides an incomplete description of the 
optomechanical scheme. 

Let us consider a simplified system with a single me- 
chanical mode, representing the movable mirror (with 



mass m and frequency ojm) of a coherently driven optical 
cavity. The optomechanical coupling between the mirror 
and the cavity field is realized by the radiation pressure. 
In the adiabatic limit in which the mirror frequency is 
much smaller than the cavity free spectral range c/2L 
{L is the cavity length) one can focus on one cav- 
ity mode only (with annihilation operator 6, frequency 
ujc and cavity decay rate 7c) because photon scattering 
into other modes can be neglected. This adiabatic regime 
implies ujm and therefore the generation of pho- 

tons due to the Casimir effect, retardation and Doppler 
effects are completely negligible. The cavity mode is 
driven by a laser field with input power p and frequency 
luq i-^c- The dynamics of the system can be described 
by the following set of coupled quantum Langevin equa- 
tions (QLEs) (in the interaction picture with respect to 
tkuob^b) 

Q{t)=uj,nP{t), (la) 

P{t) = ^LUrnQit) - lrnP{t) + Gh'^ {t)h{t) 

+ W{t)+f{t), (lb) 



iuc + icjo - y ) b{t) + 2iGQ{t)h(t) 
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where Q and P are the dimensionless position and mo- 
mentum operator of the movable mirror, 7™ is the me- 
chanical damping rate, G — {uJc/ L)yjh/ 2muj,m is the 
coupling constant, f{t) is the classical force to be de- 
tected, and E = \J p^d'hooo- The noise terms in the 
QLEs are given by the usual input noise operator bin(t) 
|lO| , associated with the vacuum fiuctuations of the con- 
tinuum of electromagnetic modes outside the cavity, and 
by the random force W(t) describing the Brownian mo- 
tion of the mirror caused by the coupling with other 
internal and external modes at the equilibrium tem- 
perature T. The optical input noise correlation func- 
tion is {hn{t)bl„{t')) = 6{t - t') 0, while that of 
the quantum Langevin force W{t) is given by [ pl]jT2[ 

{W{t)W{t')) =^ [-ir,/2T,L0m)[Tr[t-t')-iT,{t-7i 

where Jv(t) ~ dkjj lo cos{u;t) coi\\ {hu /2k bT), Ti{t) — 
dujuj sin{ujt) with ru the frequency cutoff of the reser- 
voir spectrum. The QLEs (1), supplemented with the 
above correlation functions, provide an exact description 
of the system dynamics, valid at all temperatures Q. 
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In standard interferometric applications, the driving 
field is very intense. Under this condition the system 
is characterized by a semiclassical steady state with the 
internal cavity mode in a coherent state and a 

new equilibrium position for the mirror, displaced by 
G|/3p/a;„i. Then the dynamics is well described by lin- 
earizing the QLEs (1) around the steady state, and if we 
rename with Q{t) and b{t) the operators describing the 
quantum fluctuations around the classical steady state, 
one gets 

Q{t) = UJraP{t) , (2a) 

P(i) = -u^Qit) - 7„,P(t) + Gf3 [b{t) + 6t(t)] 

+ Wit) + fit), (2b) 

Ht) = (-y - Kt) + 2iGl3Qit) + ^Anit) , (2c) 

where we have chosen the phase of the cavity mode field 
so that /3 is real and A — LOc — LOo — G'^0^/^m is the cavity 
mode detuning. We shall consider from now on A = 0, 
which can always be achieved by appropriately adjusting 
ujQ. In this case the dynamics becomes simpler, because 
only the phase quadrature Yit) — i \b\t) — 6(i)] /2 is af- 
fected by the mirror position fluctuations Qit), while the 
amplitude field quadrature Xit) = [6(t) + b\t)]/2 is not. 
The large cavity bandwidth limit 7c 3> Gj3 is commonly 
considered, and in this case the cavity mode dynamics 
adiabatically follows that of the mirror position and it 
can be eliminated, i.e., 
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Qit) 
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where Yin it) 



In this limit, the mov- 



blit)~hnit)_ 

able mirror can be used as a pondcromotive meter to de- 
tect a weak force fit) acting on it [|l3|, which will be pro- 
portional to the displacement from the equilibrium posi- 
tion Qit). The measured quantity is the output homo- 
dyne photocurrent, Youtit) = '^.^cVYit)- ^Y^t) |o), 
where -q is the detection efficiency, and Y^it) is a gener- 
alized input noise, coinciding with the input noise Yinit) 
in the case of perfect detection rj — 1, and taking into ac- 
count the additional noise due to the inefficient detection 
in the general case rj < 1 This generalized input 

noise can be written as Y^^it) ^ i \b\^it) — &,,(t)], with 
(brfit)h\^(t')) = 6it—t'), and it is correlated with the input 
noise &i„(i) according to {bin it) bl^ it')) = (&,,(t)fe|„(t')) = 

VvSit-t') 

The output of the homodyne measurement may be 
used to devise a phase-sensitive feedback loop to control 
the dynamics of the mirror. The effect of the feedback 
loop has been described using quantum trajectory the- 
ory and the master equation formalism in Ref. [Q, 
and a classical description neglecting all quantum fluc- 
tuations in Ref. j^. Here we shall use a more general 
description of feedback based on QLEs for Heisenberg 
operators, first developed in Ref. |lq] and generalized to 



the non- ideal detection case in Ref. In particular 

we shall give the first quantum description of the cold 
damping scheme 0. The adoption of a fully quantum 
treatment is justified by the fact that, as we shall see, 
in the presence of feedback the radiation quantum noise 
has important effects, especially at low temperatures. 

In the proposal of Ref. Q| , feedback induces position 
shifts controlled by the output homodyne photocurrent 
Youtit)- This is described by an additional term in the 
QLE for a generic operator Oit) given by EJ] 



/7c 



Ofbit)^t^Youtit 



r) [gs,Pit),Oit)] 



(4) 



where r is the feedback loop delay time, and gsc is the 
feedback gain. The feedback delay time is always much 
smaller than the typical timescale of the mirror dynamics 
and it can be neglected. After the adiabatic elimination 
of the cavity mode, the mirror QLEs become 



Qit) = iu,nPit) + 4G/35,,Q(t) 



Ym + 9scVl^^nit), 



Pit) 
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-LOmQit) - 7mP(t) + —=X^nit) 



Wit) + fit), 



(5a) 
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where X,„(i) = 6|„(t) + 6,„(t). 

In the cold damping scheme of Ref. |^ feedback is pro- 
vided by the radiation pressure of another laser beam 
intensity-modulated by the time derivative of the out- 
put homodyne photocurrent, and therefore one has the 
additional term in the mirror QLEs 



Ofbit) 



i dYoutit-r) 



dt 



[gMt).Oit)] 



(6) 



Again, in the zero delay limit t = 0, and adiabatically 
eliminating the cavity mode, one has 



Qit) = LOmPit) 



(7a) 



2G/3. 



Pit) = -LOn^Qit) - 7„P(t) + -^X.nit) + Wit) + fit) 



^Gpgcd 

Ic 



Qit) 



gcd 



Y^nit) 



^l(t). (7b) 



We have introduced the new input noise Yi„it), with 
the correlation function {Yinit)Yinit')) = —Sit — t'), and 
the same holds for the associated generalized input noise 

The two sets of equations (g) and (|7|) show that 
the two feedback schemes are not exactly equivalent. 
However it is possible to see that they have very 
similar physical effects on the mirror dynamics con- 
sidering, for example, the Fourier transform of the 
mechanical susceptibility in the two cases, that is, 
Xsc(w) = Wm/ [t^m + 51 7m -uj'^ + (7m + Qi)] in the 
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"stochastic cooling" feedback scheme of Ref. |Q] (51 = 
-AGPgsc), and Xcrf(w) = ojm/ [w^ - uj'^ + iuj (7™ + 52)] 
in the cold damping feedback scheme of Ref. |^ (52 = 
AG(3gcd^mllc)- These expressions show that in both 
schemes the main effect of feedback is the modification 
of mechanical damping 7^ — )■ 7„i + gi. In the stochastic 
cooling scheme one has also a frequency renormalization 
uf^ uf^ + (?i7m, which is however usually negligible 
since 7™ <C iOm- If the gains gi are appropriately chosen, 
one has a significant increase of damping and this in- 
crease is at the basis of the cooling mechanism proposed 
in Ref. [|j and realized in Ref. ||]. In fact, the mechani- 
cal susceptibility at resonance is inversely proportional to 
the damping coefficient and, in the presence of feedback, 
the oscillator becomes much less sensitive to the thermal 
noise, yielding a complete suppression of the resonance 
peak in the noise power spectrum. 

The classical treatment of Ref. ||^ (see also is 
equivalent to replace the QLEs (||) and (0) with classi- 
cal stochastic equations in which the back action noise 
2G(3Xin{t)/ and the feedback- induced noise terms 
proportional to Fi„(t), Y^^{t), Yi„(t) and Y^Ht) are ne- 
glected. As a consequence, within the classical approach, 
the only effect of feedback is damping renormalization 
7m ^ 7m + gi, while the increase of noise for increas- 
ing feedback gain, due to the presence of the feedback- 
induced terms, is completely missed. This gives the 
wrong impression that, at least in principle, unlimited 
cooling can be achieved for increasing feedback gain. 
With a quantum treatment, the tradeoff between damp- 
ing renormalization and feedback-induced noise leads to 
the existence of an optimal feedback gain, corresponding 
to the best achievable cooling. This limit on thermal 
noise suppression would be present even in the case of 
an ideal feedback loop with no electronic noise. It is a 
manifestation of quantum effects, showing the conceptual 
difference with a purely classical description of cooling. 

When the classical force we want to detect is char- 
acterized by a characteristic frequency, say, f{t) — 
fo exp [—{t — ti)^/2cr^] cos (i^ft), spectral measurements 
are commonly performed and the detected signal is 



S{lu) = 



+ 00 



dte 



(8) 



where Ft^ (t) is a "filter" function, approximately equal 
to one in the time interval [0, Tm] in which the spectral 
measurement is performed, and equal to zero otherwise. 
Using Eq. (||) and the input-output relation, the signal 
can be rewritten as 



8Gl3f] 



-t-00 



duj'x{uj')f{uj')FT^X^ - uj') 



, (9) 



where f{uo) and Ft^ (w) are the Fourier transforms of the 
force and of the filter function, respectively, and x('^) 
is equal to Xsc(w) or Xcdi'^), according to the feedback 
scheme considered. The noise associated with the mea- 
surement of the signal of Eq. (H) is given by 



r p+00 r+oo 

N{u;) = l dtFr^it) dt'FT^{t')t 

K.J —00 J —00 

^{Yout{t)Yout{t'))f=oY'\ 



j(t-t') 



(10) 



where the subscript f — means evaluation in the ab- 
sence of the external force. Using again Eq. (|^) and the 
input noises correlation functions, the spectral noise can 
be rewritten as 



iV(.).|i^ 
I 7c 



+00 



+ CXD 



xe 



C{t,t') 



dtFr^ {t) 
dtFr^ {tf 



dt'Fr^it') 
1/2 



(11) 



where C{t,t') = {Q{t)Q{t') + Q{t')Q{t))/2 is the sym- 
metrized correlation function of the oscillator position. 
Spectral measurements are usually performed in the sta- 
tionary case, that is, using a measurement time Tm much 
larger than the oscillator relaxation time, Tm ^ l/7m- 
In this limit one has Ft^ (i) ~ 1, Vt, and the signal 5'(aj) 

simply becomes S{bj) — 8G/3ri x(ijj)/(tj) jl-K^f^c- The 

oscillator in this case is relaxed to equilibrium and G(t, t') 
in Eq. ( pT] ) is replaced by the stationary correlation func- 
tion C{t — t'). Defining the measurement time T,,, so that 
Tm — J dtFT^(tY , Eq. ( |rT| ) assumes the usual form 



N(uj) 



Ic 



1/2 



(12) 



where Nq{lu) ~ J dte^'^'^'^ C (t) . In this stationary case, 
the SNR can be calculated for both feedback schemes, 
and one gets 



2o2 




for the cold damping scheme, and a similar expression 
for the stochastic cooling scheme ((j^w^ is replaced by 
gf{uj^ + 7^) and Xcd by Xsc)- In both cases, feedback 
does not improve this stationary SNR at any frequency, 
due to the gf term. This is not surprising because the 
effect of feedback is to decrease the mechanical suscepti- 
bility at resonance, so that the oscillator is less sensitive 
not only to the noise but also to the signal. 

However, in the case of an impulsive force with a time 
duration a ^ l/7m, the force spectrum could still be 
well reproduced even if a much smaller time Tm, such 
that cr <ti Tm ^ l/7m, is used. This corresponds to a 
nonstationary situation because the system is far from 
equilibrium during the whole measurement. In this case, 
the noise spectrum is very different from the stationary 
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form of Eq. (|1J) and it is mostly determined by the ini- 
tial state of the oscillator. It is therefore quite natural to 
devise a strategy in which the feedback cooling scheme is 
applied before the measurement, so that the state at the 
beginning of the measurement is just the cooled, equilib- 
rium state in the presence of feedback, and turn off the 
feedback during the spectral measurement. In this way 
the noise remains small during the whole measurement 
because the heating time l/7m is much larger than Tm, 
while at the same time the signal is not significantly sup- 
pressed because the mechanical susceptibility is just that 
in the absence of feedback. One expects that as long as 
the measurement time is sufficiently small, Tm ^ l/7m, 
the SNR for the detection of the impulsive force (which 
has now to be evaluated using the most general expres- 
sions (H) and (11)) can be significantly increased by the 
above nonstationary strategy. 

This scheme can be straightforwardly applied when- 
ever the "arrival time" ti of the impulsive force is known: 
feedback has to be turned off just before the arrival of the 
force. However, the scheme can be easily adapted to the 
case of an impulsive force with an unknown arrival time, 
as it is the case of a gravitational wave passing through 
an interferometer. In this case it is convenient to re- 
peat the process many times, i.e., subject the oscillator 
to cooling- heating cycles. In fact, cyclic cooling has been 
proposed, in a qualitative way, to cool the violin modes 
of a gravitational waves interferometer in ||l^ . Feedback 
is turned off for a time T„i during which the spectral 
measurement is performed and the oscillator starts heat- 
ing up. Then feedback is turned on and the oscillator is 
cooled, and then the process is iterated. Cyclic cooling 
is efficient if the cooling time Tcooh which is of the order 
of (7m + gi)~^-, is much smaller than Tm- This is verified 
at sufficiently large gains and it has been experimentally 
proved in iQ. In the impulsive force limit a ^ T,„, the 
performance of the scheme is well characterized by a time 
averaged SNR, i.e., 

^(-)^J^ij- Tdt,'-^. (14) 

Jm + -I cool Jo JV (UJ) 

This average SNR can be significantly improved by cyclic 
cooling, as it is shown in Fig. 1, where SNR{iL!) is plot- 
ted (full line) in the case of the impulsive gaussian force 
chosen above: one has an improvement by a factor 10 at 
resonance with respect to the no feedback case (dashed 
line). Parameter values are those of Ref. except that 
we have considered T — 4K, and the optimal value of 
the feedback gain 52 — 3057^ maximizing the SNR. The 
dotted line refers to the classical SNR{ui) calculated ne- 
glecting all quantum noises: the noise is underestimated, 
with a 15% error in the SNR at resonance. This shows 
that quantum noise has an appreciable effect already at 
liquid He temperatures and that a fully quantum treat- 
ment is needed for a faithful description of the physics. 
It is possible to see that at T = 300K, at the correspond- 
ing optimal feedback gain, there is instead no appreciable 
difference between the classical and quantum predictions, 



and that SNR{u!m) becomes 16 times larger than that 
with no feedback. 
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FIG. 1. The averaged spectral SNR of Eq. ( |l4| ) in the pres- 
ence of the cold damping feedback scheme with 32 ~ 82.4 
kHz in the full quantum treatment (full line), in the classical 
approximation (dotted line), and without feedback (dashed 
line). Other parameter values are: u>f — ujrn = 11-7 MHz, 
7™ = 270 Hz, G^P^f-yc = 10.3 Hz, T = 4 K, ry = 0.99, 
a = 3.7 /isec, Tm = 11 /xsec, Tcooi = 2 ^sec. 



In conclusion, we have presented the first quant_um de- 
scription of the cold damping feedback scheme 
we have shown how the cooling schemes of Refs. 
be used, within an appropriate nonstationary strategy, to 
improve the detection of weak impulsive forces. 
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